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(54) Digital moire profilometry 

(57) Projectors (e.g. 24) and cameras (e.g. 20) 
disposed about an object under test obtain digital 
moire information of overlapping aspects (views) of 
the object The system is calibrated using a 
calibration pattern (10) of straight lines on a flat 
reference surface and a projected pattern (26) of 
lines which are perpendicular to the lines of the 
calibration pattern and contain marker fringes 
between predetermined lines. Each aspect is 
calibrated individually by stepping the calibration 
pattern along the axis of each camera so as to 
obtain a plurality of tables which relate phase 
information to distance in each of a plurality of 
planes spaced successively closer to the camera. 
The calibration tables define a calibrated volume in 
space between the cameras and projectors. When 
the object under test is located in this volume, the 
video signals from the camera are converted into 
digital moire information. This information is 
directly related to the calibration information and is 
converted with the aid of the calibration information 
into absolute distance measurements along the 
perimeters of cross-sections through the object. 
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.method and apparatus for digital moire 
p ro f i lom £7 r y calibrated for accurate 
conversion of phase information into 
dis:a:.'j:; ^asur^msnis in a' plurality of directions 



5 The present invention relates to digital moire 

interf eromet ry , and particularly to a system and method 
using a calibrated .volume to provide absolute 
measurements of distance from digital moire phase 
information ootained from an object under test, 
10 regardless of phase ambiguities in such information. 

The invention is especially suitable for use in 
machine vision systems for measuring the contour of 
three dimension objects through selected cross-sections 
thereof. Such measurements may be used for inspection 
15 of objects as .veil as for defining their location in- 
space . 

Di^n. ^1 moire techniques are known where a line 
grating is projected on an object under test. Digitized 
vivieo information is obtained with a camera viewing the 

20 oojec:. undc*r L^st- at an angle with respect to the axis 
of the optics which project the grating. Phase 
information is derived from the digitized video using a- 
spatial synchronous algorithm (by convolving the 
digitized video information from the camera with 

25 orthogonal reference sine and cosine functions). See 
K.H. Womack, " I nt er f eromet r ic Phase Measurement Using 
Spatial Synchronous Detection", Optical Engineering, 
Vol. 23 No. 4, Page 391, July/August 1984. Phase 
information obtained as a result of such measurements is 

30 only relevant with respect to a reference plane and 
becomes ambiguous at distances greater than a half 
period of the projected grating as it is viewed at the 



BNSDOCID: <GB_2204397A^_I_> 




8NSDOCID: <GB_2204397A_I_> 



22043S7 




BNSDOCID: <GB 2204397 A 1 > 




BNSOOCID: <GB__2204397A W _I_> 




BNSOOCID: <GB 220*397 A_l _> 




BNSDOCID: <GB 2204397A_I_> 



\3\\% 

22043 




Camera Video #1 Camera Video #2 Camera Video #3 



Fig. 12A 



A 




X(Pixels) 

Fig. 12B 




X distance 

Fig. 12C 



BNSOOCID. <<3B 22Q4397A I > 



2204387 




N 



X 



g 

-a 

3 I 

life 



l 

life 
UH7 





i i 


i 4 i 




\ ST. 





Hi 
CO 



& << 

o 

>< ^ 



T 



X 



T 

N 




CO 



BNSOOCID: <GB__2204397A__L> 



22043 




BNSOOCID: <GB_2204397A__I_> 



2204307 




Fig. 10 



BNSDOCID: <GB 220 4397 A_l_> 




Fig. 9 



BNSOOCIO. <GB 2204397A I > 




BNSOOCID. <GB 2204397A_I_> 



2204397 



Convolution 
window for pixel 
(X=122,Y=Wi) 



60 

Etched lines 



! Window 

Y=Wi 




/ Projected lines 






or 

o 

3 



Fig. 7A 



Cumulative Phase 



116 



■ X • 




122 



128 



134 



Z=2 



,Z=1 



140 Pixel # 



• DATA POINTS: Shown as pixels 1 22, 1 28,134 in Fig 7A 

X STORED VALUES: Interpolated from DATA POINTS, these 
values are always stored for pixels 128 ± 6n 



Fig. 7B 



BNSDOCIO: <GB 22Q4397A I > 



611% 



2204387 




Fig. 6A 




Plane tt . 



Fig. 6B 



BNSOOCID: <GB___2204397A_I_> 



22043S7 



z. 









■ 














v 3£r£ 




.""J.**-' »■ 




* ■ 























Y=Wi 

^Convolution Window 
for Pixel 

(X-131;Y-Wi) 



Fig. 5A 



Phase /Distance 



1080°/ 0.75" - 




720° / 0.5" - 




360°/ 0.25"- 





X Pixel 



-720°/ -0.5" 
-1080° /-0.75" 

Fig. 5B 




BNSOOCID: <GB 22Q4397A I > 



f 



2204397 




BNSOOCID: <GB__2204397A_I_» 



31 \% 



22043S7 



Center Line 



!7I 



▼ 



7 



60 



Fig. 3 



BNSDOCID: <GB 2204397 A I > 



v. 



2/ IS 



22 Vision Computer 



#1 



24 




120 



X-Calibration 



Part Detect 
Pass /Fail 

Measurement 
Data 

Strobe lights 

Constant Lights 

Camera Syncs 

Camera Video 

Translation 
Rotation 



Z-Calibrauon 




V 



il_2 



Fig. 2 



Projector 




•Camera 



BNSDOCID: <GB 2204397 A_l_> 



jl<6 



2204337 



26Ruling 

24 Projector 




46 Continuous lighting 



20 CCD Camera 



Sync / Video 



Continuous 
Strobe 



-jA Projector 
Optical Axis 




Rotation Axis 



■ - 10 Calibration Target 



|l2Translation 



J16J 

4l4Rotation 




Fig. 1 



BNSOOCID: <GB I > 



POOR 
QUALITY 



UK Patent Application m GB m 2 204 397JV 

(43) Application published 9 Nov 1988 



(12) 



(21) Application No 8809852 

(22) Date of filing 26 Apr 1988 



(30) Priority data 
(31) 045110 



(32) 30 Apr 1987 (33) US 



(71) Applicant 

Eastman Kodak Company 

(incorporated In USA-New Jersey) 

Eastman Kodak Company, New Jersey, 
United States of America 

(72) Inventors 

Kenneth H Womack 
Brian J Kwarta 
David H Outterson 
James R Red a 

(74) Agent and/or Address for Service 
Swann Elt & Company 
31 Beaumont Street, Oxford, OX1 2NP 



(51) INT CL 4 

G01B 11/24 

(52) Domestic classification (Edition J): 

G1A A2 A3 C1 C3 C5 CA D10 D1 G13 G15 G16 
G17 G18 G1 G3 G7 G9 P10 P14 P17 P4 R7 S10 
S3 T14 T26 T3 T4 T7 T8 
U1S 1881 2136 2145 2148 2185 G1A 

(56) Documents cited 
None 

(58) Field of search 
G1A 

Selected US specifications from IPC sub -classes 
G01B G01N 



(54) Digital moire profilometry 

(57) Projectors (e.g. 24) and cameras (e.g. 20) 
disposed about an object under test obtain digital 
moire information of overlapping aspects (views) of 
the object. The system is calibrated using a 
calibration pattern (1 0) of straight lines on a flat 
reference surface and a projected pattern (26) of 
lines which are perpendicular to the lines of the 
calibration pattern and contain marker fringes 
between predetermined tines. Each aspect is 
calibrated individually by stepping the calibration 
pattern along the axis of each camera so as to 
obtain a plurality of tables which relate phase 
information to distance in each of a plurality of 
planes spaced successively closer to the camera. 
The calibration tables define a calibrated volume in 
space between the cameras and projectors. When 
the object under test is located in this volume, the 
video signals from the camera are converted into 
digital moire information. This information is 
directly related to the calibration information and is 
converted with the aid of the calibration information 
into absolute distance measurements along the 
perimeters of cross -sections through the object. 
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camera. Techniques for removing the ambiguity have 
involved taking into account the geometric relationship 
of the projector and the camera ( t r iangulat ion ) to 
obtain a reference plane and stepping of the grating to 
5 different positions along the path of the light 
projected on the object in order to obtain phase 
information representing the profile of objects which 
are large with respect to the grating. See U.S. Patent 
No. 4,641,972 issued February 10, 1987. 

10 Other schemes for resolving the ambiguity in 

making contour measurements have involved use of light 
of different wavelengths, synchronously moving the 
object under test and the projector and/or camera optics 
or providing other special geometric orientation of the 

15 optics. Various examples of such schemes are described 
in the following patents, 'which deal with contour 
measurements using inter feromet ry techniques, and which 
also describe complex computational processes used in 
deriving the phase information from the video signals: 

20 U.S. Patent Nos . 4 , 456 , 339 issued June 26 , 1984 ; 

4,498,770 issued February 12, 1985; 4,641,972 issued 
Feoru.iry 10, 1987 (the last three patents involve 
multiple frequency systems); 4,131,365 issued 
December 26, i973; 4,391,526 issued July 5, 1933; 

25 4,457,625 issued July 3, 1984; and 4,627,734 issued 

December 9, 1986 (the last four patents involve moving 
optics/object systems); 4,009,965 issued March 1, 1977; 
4,139,304 issued February 13, 1979; 4,340,306 issued 
July 20, 1982; and 4,387,994 issued June 14, 1983 (the 

30 last four patents involve contour measurements using 
spatial orientations of the projection and detection 
opt ics ) . 

In the various systems and techniques which 
have heretofore oeen proposed, measurement accuracy is' 
35 limited oy the alignment of the various optical elements 
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{zhe projectors and the cameras or other photoelectric 
detectors), aberrations in the lenses and the proximity 
of the reference plane to the surface of the object 
under test from which the projected light is reflected. 
In addition, the measurements provided by such systems 
are only the relative phase with respect to the 
reference plane. 

It has been discovered in accordance with the 
invention that an absolute measurement of surface 
coordinates (dimensions in a plurality of coordinates) . 
can be obtained by calibrating the phase information 
with respect to a calibration grating (hereinafter also; 
referred to as a ruling) in a plurality of planes at 
different incremental distances from the detector (e.g. 
a video camera). For surface coordinates of a solid 
oo^ect successive aspects (which may also be considered 
as views) of the calioration grating are. used. The 
information ootainsd from the calibration is stored in 
tables to provide a calibrated volume, such that 
measured phase values can be converted directly into 
distance measurements (to dimensional coordinates in a 
reference coordinate system). Since the calibration 
takes into account the geometric orientation (alignment) 
of the optical elements and the aberrations and other 
distortions of the optical elements (of the projector 
and video camera) and any distortion in the sensors (the 
video camera), the absolute measurements are 
automatically compensated for any such misalignments, 
aberrations and other distortions. 

A profile of a three dimensional object can be 
provided by utilizing a plurality of camera and 
projector arrangements which cover different aspects of 
the object. The volume in space which is covered by 



- 4 - 



each camera and projector arrangement is calibrated so 
as to provide a three dimensional volume in« which the 
object under test can be located. The measurements of 
the distances defining the profile of the object under 

5 test through a cross-section thereof is obtained by 
combining the information obtained from each of the 
several aspects which are covered by each 
camera/projector arrangement. 

It is a feature of the invention that errors in 

10 the precision of the measurements are determined by the 
accuracy of the calibration grating and the mechanical 
stage on which it is located. These errors are 
therefore controllable to the degree of precision 
required, for example a fraction of a thousandth of an 

15 inch. 

Accordingly, it is the principal object of the 
present invention to provide an improved system of and 
method for prof i lome t r y , both of surfaces and three 
dimensional objects, which provides absolute 
20 measurements of the dimensions thereof with precision 
and which is not limited by misalignments, aberrations 
or distortions in the optical elements or sensors of the 
system. 

It is another object of the present invention 
25 to provide an improved system for making profile 

measurements both of surfaces (two dimensional objects) 
and three dimensional objects with precision and without 
the need for scanning or shifting the position of the 
detectors with respect to the object under test to 
30 provide measurements of large objects. 

It is a still further object of the present 
invention to provide an improved method and system for 
precision measurements by i nter f eromet ry which 
eliminates ambiguities in phase measurements and 
35 provides absolute measurements of distance defining the 
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contour of objects under test, whether in two or three 
dimensions . 

It is a still further object of the present 

invention to provide an improved contour measurement 
5 system which utilizes spatial phase measurements to 

correct for the effect of temperature fluctuations on 

the geometry of the system. 

It is a still further object of the present 

invention to provide an improved method of and system 
10 for profilometry in which a ruling is projected upon an 

object under test wherein the ruling provides a fringe 

pattern encoded, as with marker fringes (marked ruling 

lines) between predetermined lines of the ruling, so as 

to provide information as to absolute distances of the 
15 surface under test from a reference position 

irrespective of ambiguities in the phase information. 

The foregoing and other object features and 

advantages of the invention 

are illustrated by and 
20 apparent from a reading of the following description in 

connection with the accompanying drawings, listed below 

and given by way of example: 

FIG. 1 is an elevational view schematically 
showing a system embodying the invention. 
25 FIG. 2 is a plan view schematically 

illustrating a system embodying the invention which is 
capable of measuring the profile of an object in three 
dimensions. 

FIG. 3 is an elevational view of the 
30 calibration target which is used in calibrating the 
system illustrated in FIGS. 1 and 2. 

FIG. 4 is a view of a variable frequency Ronchi 
grating used in the system shown in FIGS. 1 and 2; the 
grating having marked lines which provide marker fringes 
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between predetermined rulings (the ruling along the 
horizontal center line of the grating and the next 
vertically displaced ruling and rulings eight periods 
below tne marked ruling) so as to enable calibration of 
5 a large volume where the object under test might be 
located. 

FIG. 5A shows a portion of a horizontal window 
at a height W. ; the window containing a pattern of 
lines of the calibration grating (FIG. 3) about the 
10 center line thereof as imaged on the sensors of the 

camera. The horizontal window reflects a row of pixels 
at y=H i# which is calibrated. There are 32 
calibration windows. The figure is calibrated in 
accordance with the location of the pixels sensed along 

15 a horizontal line — 256 pixels with pixel number 128 at 
the center line and 127, 128 pixels on either side of 
the centerline. 

FIG. 5B is a curve illustrating the variation 
in the phase measured at the pixel locations on the 

20 right side and left side of the center line, 

respectively. The correspondence between distance and 
phase being shown on the ordinate of the curves. 

FIG. 6A is a diagram showing the window 
illustrated in FIG. 5A with the projected Ronchi grating 

25 superimposed thereon and showing the location of the 

marker fringe along the center line of the Ronchi 

grating when the calibration target (FIG. 3) is stepped 

to a certain location, Z closer to the camera than an 

n 

initial location of the calibration target at Z Q . 
30 FIG. 6B is a curve showing the cumulative phase 

(the phase measurement at the window to the phase 
measurement at the marker fringe plus the number of 360° 
phase transitions therebetween) for each plane at which 
the calibration grating (FIG. 3) is located, plane Z n 
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oeing furthest from the camera and plane Z 3 being 
closest to the camera. 

FIG. 7A is a diagram similar to FIG. 6A showing 
the location of a convolution -window which is used to 
measure the phase in the horizontal direction along a 
window W... These phase measurements being used for 
calibration in each of the planes Z Q to Z in which 
the calibration grating (FIG. 3) may be located. z 
being the plane at the origin with higher numbered 
planes being successively closer to the camera. 

FIG. 7B is a series of curves showing the 
cumulative phase between measurements with the 
convolution window at different pixel positions for each 
of three successive planes Z Q , z l and Z 2 ; the 
15 remaining planes not being shown to simplify the 
illustration. 

FIG. 8 is a diagrammatic view in an elevationai 
orientation as shown in FIG. 1. fig. 8 shows the 
location of the projected Ronchi grating in successive 
20 calibration planes 2 n and Z (n+i) and also showing 
the location of an intermediate plane. The figure 
illustrates how the accumulative phase measurements made 
at the calibration planes may be interpolated to 
determine the phase measurement and the distance at an 
25 intermediate plane between the calibration planes. 

FIG. 9 is a view in a horizontal plane of the 
calibration planes Z n and Z r+1 and also another 
calibration plane Z^ further from the camera than 
2 n and showing the location of the surface of the 
30 object under test with respect to these planes so as to 
illustrate how the phase measurements at each pixel 
location in each plane through a cross-section of the 
object is mapped to the distance of the object, directly 
from the phase information using the z calibration 
35 information i ch is plotted in FIG. 7B. 
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?IG. 10 is a diagram showing the relationship 
between absolute phase (0) and unwrapped phase U{.fT) 
which takes into account successive 360° transitions. 
The curves show how absolute phase measurements are 
5 ambiguous, and how the ambiguity can be removed through 
the use of unwrapped or cumulative phase. 

FIGS . 11A and B together constitute a 
functional block diagram of the profilometer system as 
implemented with the three sets of cameras and 
10 projectors which cover successive overlapping aspects of 
an object under test and the computer system as shown in 
FIGS . 1 and 2; 

FIGS. 12A, B and C show: in FIG. 12A the 
images of a cylindrical object of elliptical 
15 cross-section imaged at each of three different cameras, 
120° apart as shown in FIG. 2; in FIG. 12B the variation 
in phase in the x direction cross a horizontal 
window W. ; and in FIG. 12C the distances in the z and 
x directions into which the phase information is 
20 converted through the use of calibration information. 

FIG. 13 is a functional block diagram 
illustrating in greater detail how the phase 
calculations ace performed in the system illustrated in 
FIGS. 11A and 1 1 B ; 
25 FIG. 14 is a functional block diagram 

illustrating the operation of the calibration system to 
obtain the tables in memory in the system shown in 
FIGS. 11A and 11B. 

FIGS. 15A, 15B and 15C are functional block 
30 diagrams illustrating in greater detail the calibration, 
marker fringe table generation, and z calibration, 
respectively. 
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Referring to FIGS. 1 and 2 there is shown, the 
profilometer system in accordance with an embodiment of 
the invention set up for calibration. A calibration 
fixture 11 has a calibration target 10 mounted on a 
5 stage 12 and 14 which permits the target to be rotated 

about a vertical axis which is indicated as the "y" axis 
and translated along the "z" axis. The stage includes 
plates 16 and 18 mounted on a table 14 which provides 
for rotation. The plate 18 is translatable in the z 
10 direction. The surface of the plate 18 is leveled with 
reference to gravity. After calibration the calibration 
fixture 11 is removed and the objects under test are 
located in the calibrated volume. The objects, for 
example, may be bodies for which the contour of a cross 
15 section is to be determined. These bodies may be moved 
through the calibrated volume on a conveyor belt which 
is also leveled to gravity so that the axis of the body 
is approximately along the y axis. 

The translation stage is designed to enable the 
20 calibration target to move through 32 planes, Z 0 _ 31 . 
The Oth plane is through the y axis and the other z 
planes are a fraction of the period of a line 

pattern from a Ronchi ruling which is projected on the 
calibration target. This pattern in an exemplary 
25 embodiment may have a period of 1/4 inch. The 

translation steps may then suitably be 1/8 inch. 

While the steps may be equally spaced, their 
spacing is not critical as long as they are less than 
half the period of the ruling as projected on the 
30 calibration target. The resolution of the calibration 
and also of the profile measurements depends upon the 
period of the ruling. It has been found that, with a 
period of approximately 1/4 inch, measurement accuracy 
to + 0.001 inch, is obtainable. 
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The calibration target is rotatable to three 
positions 120° apart as shown in FIG. 2. In each of 
these positions, the calibration target, which presents 
a flat surface on which the vertical- pattern shown in 
5 FIG . 3 appears, is located opposite a projector camera 
assembly. These assemblies are indicated as PC#1, PC#2 
and PC83 in FIG. 2. All of these assemblies are 
identical. PC#1 includes a video camera which is 
preferably a charge couple device or CCD camera which 

10 has a CCD sensor array. The array may be a matrix of 
about 320 horizontal and about 250 vertical elements. 
The video signals from these elements are digitized in a 
computer system 22 which is connected to the camera 20. 
The video image is digitized into a matrix of 256 

15 horizontal and 240 vertical elements (pixels) as 

represented by bytes of digitized information from .the 
CCD sensor elements. Each of these elements of 
digitized information corresponds to a pixel of the 
field viewed by the array. This field is focused on the 

20 array by a lens in the camera. The focus is set to the 
mean location of the surface being imaged, with a 
suitable depth of field to cover the calibrated volume. 
Since the digitized array has 256 pixels in the 
horizontal direction, the center pixel is 128, the 

25 column of 240 pixels along the center (pixel 128 in each 
row) is aligned with the y axis in the plane of the z 
axis (the zy plane). Such alignment may be accomplished 
manually using conventional positioning devices. 

The projectors each include a projector lens 

30 which focuses light projected through the variable 

period Ronchi ruling 26 on the plane of the calibration 
target. While the optical axis of the camera is aligned 
with the z axis and is horizontal, the optical axis of 
the projector is tilted at an acute angle, suitably 45° 

35 to the optical axis of the camera. The camera optical 
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axis is indicated at 28 and the projector optical axis 
is indicated at 30. These axes are aligned to be in the 
z,y plane using the positioning devices. 

The Ronchi ruling 26'is illustrated in greater 
detail in FIG. 4. It desirably has a variable period. 
When an inverting projection lens is used, the lines are 
closer at the top than at the bottom. The variable 
spacing compensates for the tilt of the projector axis 
30 so that the projected lines are parallel and equally 
spaced on the calibration target 10 when in place or 
upon objects under test. 

While the projection lens is shown 
approximately focused when the calibration target is 
slightly displaced toward the camera from the vertical 
axis at Z Q , the optimum focus is set for the location 
of the calibration target surface in the middle of its 
translation range (at step 16). The focus of the camera 
20 may be similar. The F stop of the projector and 
camera is selected to provide a depth of field over the 
entire translation range. 

The Ronchi ruling as shown in FIG. 4 contains 
at line 33, which extends along the middle of the 
grating, perpendicular to the constructive center line, 
an area which provides a marker fringe 34. Also between 
lines spaced below the center ruling 33 and intersecting 
the center line is another marker fringe area 36. The 
lower fringe 36 may be eight lines below the center 
fringe 34. Optionally, there may be additional marker 
fringes between the same lines as the lines containing 
the fringes 34 and 36 equally spaced from these lines. 
These additional marker fringe areas on opposite sides 
of the center fringe area 34 are indicated at 38 and 
40. The additional marker fringe areas on opposite 
sides pf the lower center line fringe area 36 are 
indicated at 42 and 44. The fringes adjacent to ruling 
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33 are 75% transmissive, as compared to totally 
transparent (100% transmissive), the marker fringes 36, 
42 and 44 between the lower lines are suitably 85% 
transmissive. The amount of t ransmissi vi ty and the 
relationship of the t ransmi ssi vi t ies of the fringes may 
be varied or reversed. In addition further marker 
fringes spaced downwardly or upwardly from the center 
line may be used in order to provide a larger area in 
the z direction which can be calibrated and used 
(providing a larger calibrated volume for the object 
under test). As shown in the enlargement in FIG. 4, the 
selected t r a nsmi ss i v i ty is provided by a half tone area 
where the percentage t ransmissivi ty is determined by the 
percentage of the transparent to dark areas of the. half 
tone . 

Another optional means for covering and 
caliorating over a larger area in the z direction is to 
use a reduced angle of tilt, for example less than 45° 
between the optical axes 28 and 30. The use of 
additional marker fringes for covering larger areas is 
preferred since greater resolution is obtained than 
would be the case by decreasing the tilt angle. 

Each of the camera projector assemblies PC#l-#3 
covers a separate aspect (or view) of the object under 
test and the calibrated volume; 120° aspects in this 
embodiment. Each 120° aspect is calibrated separately 
so as to cover the entire volume. FIG. 2 illustrates 
the Z Q plane for the aspect viewed by P/C#l, which is 
along the rotation axis (y) and several of the 
additional z planes each of which is successively closer 
to the camera. Exemplary z planes for the other aspects 
which are viewed by the other camera and projector 
assemblies are also* shown. 

The calibration target is shown in detail in 
FIG. 3. It has accurate "rulings which may be fabricated 



- 13 - 



by a photolithography process on a planar surface, such 
as a glass plate. The lines are suitably 1/8 inch wide 
black lines separated by a 1/8 inch white space, thereby 
providing a 1/4 inch period. 'The selection of line 
spacing is governed by the resolution required, the 
Field of View to be covered, and the receptor matrix in 
the sensor used. Since the accuracy of the calibration 
depends upon the precision of the lines, they are 
desirably made to an accuracy which exceeds the 
measurement of accuracy of the system, for example an 
order of magnitude more accurate than the desired 
measurement accuracy. The exemplary target may be a 
square plate 8 inches on a side. The flatness of the 
ruled surface of the plate is preferably an order of 
magnitude higher than the resolution of the measurements 
(e.g., in this case 0.0001 inch). It is desirably 
mounted in the calibration fixture 11 so that the lines 
are parallel to the rotation axis. The alignment of the 
camera projector assemblies is facilitated by the outer 
marker fringes so that the projected ruling pattern at 
the target will constitute lines which are perpendicular 
to the lines of the calibration target. The alignment 
is not critical since the calibration process 
compensates for errors in alignment and also for 
aberrations in the lenses and sensors of the camera and. 
projector. 

The projector is illuminated by a strobe lamp 
(not shown) which may be connected to the projector 
housing by fiber optics. A general illumination light 
source is provided on each camera projector assembly. A 
typical one of these light sources is a lamp shown at 
46. These lamps provide diffuse lighting and illuminate 
the calibration target. They may be incandescent 
lamps. The projection strobe lamps for each projector 
and the illumination lamps for each camera are 
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controlled by tne computer 22. The computer provides 
horizontal and vertical sync signals to each camera. 
Separate controi lines are provided for the projector 
strobe lights and the camera constant illumination 
5 lights 46. The computer also controls the rotation and 
translation of the stages 12 and 14 by operating stepper 
motors connected to the stages. These stepper motors 
and appropriate gearing are not shown to simplify the 
i llust rat ion . 

10 The computer processes the video information 

both for calibration and during tests on actual objects 
(run time). The computer system for both calibration 
and tun time is described hereinafter in connection with 
FIGS. 11 - 15. The calibration process utilizes the 

15 calibration target rulings and the projected Ronchi 
pattern. Their use and the operation of the system 
during calibration will become more apparent from 
FIGS. 5 - 10. 

Tne system may be compensated for temperature 

20 variations as will be explained more fully below. Also 
in order to minimize the adverse effects of environment, 
including ambient light, it is desirable that the entire 
system be enclosed in a chamber. This chamber may be 
maintained slightly above atmospheric pressure in order 

25 t0 minimize air currents and to exclude dust which might 
cause optical interference. In the event that the 
facility in which the system is used is subject to 
vibration, it is also desirable to mount the entire 
system on a table which is isolated for vibration. It 

JO is also desirable to avoid large temperature changes 
(e.g., greater than +5°F) ; smaller changes are 
accommodated by a temperature control function described 
hereinafter. The system may be enclosed in a chamber 
which is temperature controlled to prevent such large 

35 temperature changes. 
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In general, calibration is carried out to map 
the relationship of phase and distance in a plurality of 
planes (the z planes) dhich extend between the rotation 
axis and the cameras. These maps contain the phase and 
5distance relationships at a plurality of different 
heights in the y direction. Each height may be 
considered to include a window which extends in the 
horizontal direction. 

The calioration is carried out first using the 
10 cal ibrat ion target alone illuminated with the constant 
illumination from the respective lamps 46. A phase 
measurement is then made from which the distances in the 
x direction from the center of the vertical lines on the 
calibration target is obtained. Distance is correlated 
15to the positions of the pixels across the row, i.e., the 
distance of each pixel from the center pixel. Because 
the vertical lines of the calibration pattern have a 
predetermined spacing (1/4 inch in this example), the 
phase measurements are also directly correlated with 
20distance. There may be 32 horizontal windows. By 

limiting to 32 windows, the computer memory capacity is 
reduced. More or fewer windows may be used, if 
desired. The windows are located at y values 
y=136 + no, where n is an integer. y may vary from 42 
25to 228. Each window contains 32 data points, which 

again may be increased or decreased as desired. The 32 
points used are defined at pixels, x=128 +_ n 6, n is an 
integer and x is in the range from 32 to 218. Thus, 
there are 16 points (pixels) on one side of center; 15 
j0 on the other side of and at center, thereby 

defining a total of 32 data points per window. There 
are 32 calibration windows per plane, and there are 32 
calibration planes 2 Q to Z 31 . Thus there are a 
total of 32 matrices of 32 by 32 values of 
35distance/pixel information. 32 pixel positions in 32 
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windows are caliorated. X calibration, I calibration, 
and a marker fringe table in each z plane are completed 
before stepping the target to the next plane closer to 
the camera. 

5 Next the Ronchi ruling with the marker fringes 

is projected onto the caiioration target by flashing the 
strobe lamp. The relative position of the window W i 
with respect to the marker fringe is determined. Any 
other one of the horizontal windows may alternatively be 

10 used. This is done by storing phase information as to 
the location of the marker fringe and phase information 
as to the location of the window plus the number of 
lines of the projected pattern (horizontal lines) 
between the marker fringe and the window. This number 

15 of lines is determined by counting the cycles of 360° of 
phase variation between the marker fringe and the 
window. U will be appreciated that the marker fringe 
location is indicated by a dip in the pattern of 
intensity variation along the vertical. or y direction. 

20 This pattern would normally be sinusoidal because of the 
equally spaced lines of the projected grating. The 75% 
and 85% fringes can be discriminated by the difference 
in the amplitude of the intensity variations. During 
some measurements for the location of the marker fringe, 

25 the fringe 34 at the center of the projected pattern may 
move out of the field of view. This may occur at planes 
close to the camera. Then the other fringe 36 will 
appear. Since it is of different transmissivi ty , the 
detection thereof can be used in place of the central 

30 marker fringe 34 merely by adding the phase displacement 
corresponding to the number of lines of separation 
between the fringes 34 and 36. In this example there 
are eight lines or 8 times 360° of phase displacement. 

Absolute phase information between the marker 

35 fringe and the window is measured in unwrapped phase. 
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The meaning of unwrapped phase will be apparent from 
FIG. 10. The sawtooth curves indicated by the dash 
lines marked 0 represents absolute phase versus 
distance. At each period of the pattern, phase varies 
5. 360°. Unwrapped phase as indicated by the line marked 
U(,0') is the unwrapped phase. The 360° cycles of phase 
variation are cumulative, thus as the distance along the 
projected ruling between the location of the marker 
fringe and the location of the window increases, the 
10 unwrapped phase will increase. For example, consider 
the absolute phase over the distance indicated as A x 
in the figure. a x starts at 315°, passes through 
360°, and continues for another 45°. In terms of 
absolute phase measurement, a x is ambiguous because of 
15 the transition in the absolute measurement. In terms of 
unwrapped phase the ambiguity is removed. At the end of 
the ,\ x interval, the unwrapped phase measurement is 
360° plus 45- or 405°. Of course the absolute phase 
measurement is only 45°, and this measurement is 
20 ambiguous as to the location of the ^ x interval with 
respect to x. 

By measurement of the unwrapped . phase between 
the marker fringe and the window at each plane of the 32 
planes from Z Q to z 31 , a phase number corresponding 

25 exactLy to the plane is obtained. This phase number 

corresponds exactly to distance in the z direction since 
the phase number for phase Z Q , corresponds to the plane 
at the axis of rotation (the vertical axis through the 
origin - see. FIG. 1 - which is the axis about which the 

30 stage 12 rotates). The phase number for the plane Z 

corresponds to a plane 1/8 inch away from the axis in. 
the z direction, the phase number for the plane Z 
•corresponds to a plane 1/4 inch from the axis of 
rotation and so forth. A table of the phase numbers 
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therefore precisely locates all of the planes in the z 
directions. 

Maps calibrating each of these z planes in 
terms of the phase at each calibrated pixel position 
along the window are obtained by making absolute phase 
measurements at each pixel position. This phase is 
compared to the phase found for that pixel at the 
preceding plane. The difference accumulated, is 
referred to as Z-cal ibrat ion . These illuminated pixel 
positions are those pixels between the vertical lines on 
the calibration target. (Not every pixel is measured, 
because for some, the projected ruling pattern is of low 
visibility due to the dark rulings on the flat 
calibration surface of the target 10.) In general, due 
to perspective effect, the pixel at which the phase is 
determined is not the pixel at which the calibration 
values are stored. 3y a linear interpretation of 
adjacent phases, the phase at the calibration pixel may 
be found. Maps of phase measurements are made taking 
into account the phase measurements at the Z Q plane. 
At Z Q there is no preceding plane. The phase map thus 
originates there. Each calibration pixel phase (recall 
this may be interpreted from the phase of adjacent 
pixeis) begins as an offset from the center phase (phase 
at pixel 128). This compensates for aberration and 
distortion in the optics and sensors (cameras, etc.). 
The center phases are taken as zero at the Z Q plane 
and constitute the link between the marked fringe array 
and the Z-cal i bra t ion array. 

There are 32 planes of phase maps which provide 
the z calibration for each aspect, there is a total of 3 
times 32 or 96 maps for the three dimensional system 
shown in FIGS. 1 and 2. Since the x calibration 
obtained under constant illumination of the calibration 
target relates the pixel positions in each of these 
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Planes to distance, and since the z calibration maps of 
phase information permit conversion of the phase 
information to distance in an orthogonal direction 
absolute distance measurements" in an (x,z) cartesian 
system are made during run time. 

X-calibrat ion will become more apparent from 
FIG. 5A where a plurality of the lines of the 
calibration target are shown. Shown is the calibration 
procedure used to calibrate window w. of plane. Z 
There are 32 calibration windows, W., at each Z 
plane. These are the lines about the center line or 
Pixel 128 as viewed by the camera. A convolution window 
one pixel wide in the y direction and 25 pixels long in 
the x direction is moved both to the right and to the 
15 left of the center line. These phase measurements 

relating the pixels along the row encompassed by the 
window are shown in FIG. 5B in terms of unwrapped 
phase. Since each cycle of unwrapped phase (each 360° 
of" unwrapped phase), col- responds to a 1/4 inch (the 
20 distance between the centers of the white spaces between 
the lines), there is a direct correlation in the x 
direction between phase and distance. These distance 
measurements are stored in a separate region of computer 
memory for each plane, Z n - Z 31' The x calibration as 
25 ^own in figs. 5A and B is carried out under constant 
illumination and before the projector strobe is 
illuminated to project the Ronchi pattern onto the 
calibration target. 

Calibration of the window, tf. , locations 
30 using the marker fringes will be more apparent by 

referring to FIG. 6A which shows the portion of the 
window and also the portion of the pattern projected on 
the camera having the marker fringe 34. The total 
cumulative phase from the position of the marker fringe 
55 y m t to the position of the window y w . is J& 
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Total cumulative phase varies in each plane, Z . This 

n 

cumulative phase takes into account the number of lines 

of the projected pattern between the marker fringe and 

the window. The phase measurements .are made using a 

5 convolution window 1 pixel wide in the x direction and 

25 pixels wide in the y direction. The spacing between 

lines will encompass at least 3 pixels. Thus, as the 

convolution window moves upwardlv from y e to the 

mf 

position of the window y wi , there will be a sinusoidal 
10 variation in phase, with each period of such variation 

corresponding to movement of the window a distance equal 
to the period between lines of the projected pattern. 
By counting the number of these lines and adding 2 
times the numoer of lines to the phase difference 
15 between the absolute phase measurements at the position 
of the marker beam y mf and at the position of the 
window y w ., the cumulative phase is determined. This 
cumulative phase is stored for each plane in the z 
direction (Z Q to Z 31 > and correlates the distance of 
20 these planes from the rotation axis to the phase 
measurements. 

FIG. 8 shows how the position of the marker 
fringe with respect to the window varies in the 

different planes. Two planes are shown Z and Z 
„ n n+1' 

25 * n+i being the plane adjacent to which is closer 

to the camera. The window W is tho location of a 

certain row of pixels of the sensor array. Due to the 

convolution nature of this window, it must be located at 

a distance greater than 12 rows either from the top or 

30 bottom of the array (from y=12 to y=228)/ The reason 

for the appearance of 12 rather than 13 results from the 
final value of the convolution array being zero rather 
than a significant number. Because of the tilt between 
the optical axis of the projector and optical axis of 

35 the camera, the marker fringe traverses upwards as the 
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calibration plane moves closer to the camera. A 

consequence of parallax causes the window height, W. , 

to decrease as the calibration plane moves closer to the 

camera . * 

While the marker fringe 34 is shown, it will be 

appreciated that in some cases, particularly in planes 

close to the camera, the lower marker fringe 36 will be 

used and the unwrapped phase corresponding to the 8 

cycles (8 lines) between these marker fringes 34 and 36 

will be taken into account in the stored value of Of ^. 

mf 

FIG. 63 shows the variation in cumulative phase 
'° m f zn with the planes as the planes move closer to the 
camera. The plane at the- origin Z^ has a certain 
cumulative phase. . This phase continues and then 
15 reverses as the location of the marker fringe moves to 
the opposite side of the window W. . 

In FIG. 8 there is shown also between the Z 

n 

and 2 n+1 planes, a plane indicating the location of 

the surface of an object or part under test as would be 

20 the case during the run time after calibration. The 

plane of FIG. 8 is taken through the center line (a 

plane- including the axis of rotation and the optical 

axes) so that pixel 128 will be intersected by this 

plane. This pixel is indicated at p 10Q in FIG. 8. 

2S The measurement of cumulative phase Or ^ is 

mt p 

proportional to the distance from the origin of the Z 

n 

and Z^ +i planes. The phase measurements correlate to 
this distance in accordance with the following 
relationships which define the depth between the Z 

n 

30 and 2^ +1 planes (%Z) 

* Z = (0 'm£z " ^fp ,/(a mf 2ri -^fz n + l' (1 > 

where <0' m £ is the cumulative phase from the window 

W. to the marked fringe for plane Z , /f f is 

n mEz^ 
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the cumulative phase from the window w. to the marked ' 
fringe for plane 2^+1, and # mfp is the cumulative 
phase from W. at the surface of the part under test to 
the markeJ fringe. 

The distance of the surface of the object Z at 
pixel 128 is then obtained by the following expression 

Z a = Z n + %z C -a Z). (2) 

In this example the size of the step ( L Z) is 1/8 
inch. Z n is the size of the step times the number of 
the plane from the origin. For example plane 8 (Z = 
8) is I" from the origin. Accordingly, the cumulative 
phase information directly provides an absolute 
measurement of the distance from the origin of the 
surface of the part. 

Z-calibration will become more apparent from 
FIGS. 7A and 7B. The calibration is shown taking place 
at a window y . it will be appreciated that a 

VW JL 

plurality of such windows, for example 32 windows, may 
be used and located anywhere between y=12 and y=228. In 
the present embodiment, as discussed above, y=138 + n6, 
for y = 42 to 228. The convolution nature of the window, 
limiting the extent to which the calibration windows 
approach the ooarders was also discussed earlier. 

In FIG. 7A only the pattern of the calibration 
target (vertical bars) and the projected pattern 
(horizontal bars) in the vicinity of the center line 
(pixel 128) is shown. The calibration bars are marked 
"etched lines" and the lines of the projected pattern 
are marked "projected lines" in FIG. 7A. Initially the 
Z Q plane is calibrated. This is accomplished by 
moving the convolution window and making phase 
measurements at the pixel positions which are exposed 
between the etched lines. In FIG. 7A these are pixel 
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positions 116, 122, 128, 134 and 140 for piane Z^. It 
is noted that these pixel locations will vary with plane 
number, as consequence of parallax. Absolute phase 
measurements are made by the synchronous detection 
5 process which is discussed in the above-identified 

article by Womack and will be discussed further below in 
connection with the computer operations. The phase 
measurements are cumulative phase with respect to the 
Z Q plane and the phase measurements in the Z Q plane 
10 are taken into account in the cumulative phase 

measurements for all of the other 31 z planes in each 
120° aspect. 3y taking into account these phase 
measurements, misalignments and aberrations are 
automatically compensated. 
15 The phase information for the Z Q plane may be 

stored in terms of the difference in the phase 
measurements at the pixel positions outward from the 
pixel 128 along the center line to the phase at pixel 
128. Thus if the system were perfect without 
20 aberrations or misalignments, the phase measurements at 
each of these pixels would be 0. However due to 
aberrations anJ misalignments the phase measurements can 
vary. Therefore the value stored for phase at each 
pixel position will vary to accommodate these 
25 misalignments and aberrations. Recall that unwrapped 
phase was discussed in connection with FIG. 10 above- 
As the measurements of phase are made in each of the 
other planes, for example plane Z^, differences 
between the phase at each pixel location measured in 
30 plane Z^ from the measurements at corresponding pixel 
locations at plane Z Q are derived. To these 
difference values there is added accumulated phase from 
the previous plane. In the case of plane this will 
be the stored values for plane Z Q . In case of plane 
35 Z 7 the calculated values for plane Z, are then added 
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to the absolute cumulative phase values at corresponding 
pixel positions. It will be seen that the numbers 
representing the phase values increase as the planes 
move closer to the camera. In order to increase 
5 resolution these numbers may be multiplied by a scaling 
factor, for example 8k (8096) (for each 180°). 

FIG. 7B illustrates how the cumulative phase 
values changes as the planes come closer to the camera 
in the pixel positions 116, 122, 128, 134 and 140 about 

10 the center of the sensor array. The curvatures take 
into account lens and sensor aberrations and 
characteristics, parallax errors and alignment errors. 

The 2-calibrat ion phase values are used 
directly by interpolation to locate the surface of an 

15 object as shown in FIG. 9 wherein the sur'face of a part 

is depicted passing through plane Z . , Z and 

n— l n 

2 n+1 « The absolute distance calculation for each of 

these pixels can be calculated from the cumulative 

(unwrapped) phase difference in accordance with the 

20 following relationships, wherein x is the pixel number 

where the Z distances is to be measured, & is the 

phase at the pixel x, , is the phase at the x-1 

x x 

pixel adjacent to the x pixel, P is the cumulative phase 
to plane (just before — further from the camera 

25 than the unknown distance), P* is the cumulative 

phase to plane (just after the unknown distance), 

P^_ 1 is the cumulative phase in the x direction at the 
x-1 pixel, P is the cumulative phase in the x 
direction at the x pixel, Z and Z' are the distances 

jO from the origin of the calibration planes which straddle 
the distance to be measured, and U is the unwrapped or 
cumulative phase. 
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o ( ^ o) = v (jr x - ^ x _ 1 ) ( 3 ) 

P x = ? x _ j_ + U ( A ^r) (4) 

%p x = (P " p xV (p " P ' > (5) 

Z x = Z + %p x (Z - Z» ) (6) 

5 By way of an example for the surface measured at pixel 
129 (P 12 9 J the absolute distance at Z 12 $ mav De 
determined using the following calculations. 

u .o) - u t.cr 129 - ^ i28 ) ( 7 ) 

P 129 = P 128 + " ( ^ .0) ( 8 ) 

10 %p 129 - CP- p 129 )/p - pf > (9) 

z = z + %p d (z - Z' ) (10) 



As the calibration values exist for x=128 + n6 
for those pixels which are not calibrated, values are 
interpolated from those which are. For example in 
15 equations (7) to (10) the values for P. and P 1 

representing pixel 129 are interpolated from the 
calibrated values at pixels 128 dnd 134. 

RUN TIME 

Referring to FIG. 11A the system is shown as 
20 configured during run time. The computer operates by 
strobing each of the projectors in each of the three 
camera projector assemblies in succession. The camera 
input is then digitized in a 256 by 240 byte array. 
Each byte has 8 bits to indicate the gray level 
25 (variations in intensity). The control function of 
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strobing and digitization of the CCD arrays is indicated 
by the block 100. The computer stores the arrays from 
each of the 3 cameras in an area of 512 bytes by 512 
bytes in a video memory section of the computer memory. 
5 This storage may be in random access memory (RAM) in the 
computer memory. The storage arrays for each camera 
(1-3) is indicated by block 102. 

In the case of a tubular object of elliptical 
cross-section where camera #1 views in a direction 

10 perpendicular to the major axis of the ellipse, the 

video camera sees a pattern such as indicated under #1 
in FIG. 12A. Cameras #2 and #3 which are 120° apart 
have imaged thereon patterns shown under #2 and #3 in 
FIG. 12A. The calibration window which was referred to 

15 in FIG3. 5A, 6A and 7A is indicated at by way of 

example. The obtaining of the profile of the object at 
a cross-section perpendicular to the axis of the object 
which includes the window W. will be explained. It 
will be appreciated that this window may be one of 

20 several, if the entire surface profile of the object is 
to oe derived. 

Each camera utilizes a similar phase 
calculation and distance conversion subsystem 104, 106 
and 10S. The subsystem 104 for camera #1 is shown in 

25 detail. The marked fringe is first located using the 
bytes stored in the array for camera #1 in the video 
memory 102. The bytes are scanned along the column 
containing pixel 128 and the variations in gray level 
(intensity) vary sinusoidally except in the area of the 

30 marked fringe. The row containing the marked fringe 
location is located and flagged. 

The fringe is removed by scaling the gray level 
of the byte representing the marked fringe so that the 
oytes alonj pixel i28 will vary sinusoidally in 

35 amplitude. The marked fringe is removed; however, it 
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need not be removed, if it is not located in the window 
which includes the cross-section at which profile 
measurements are being made. If any of the other marked 
fringes are in any of the cross-sections, their gray 
5 level is scaled to restore the variation intensity and 
eliminate the effect of the marked fringe. 

In the next part of the phase calculation and 
distance conversion, the phase at the marked fringe and 
the phase of the pattern at the window a re determined . 
10 FIG. 6A and 8 above explained in general how this phase 
was measured and used to determine the phase at the 
object and convert that phase measurement to distance 
using equations (1) and (2). 

At this point spatial synchronous detection is 
15 used to measure the phase from the video information 
which is stored in the video memory. FIG. 13 
illustrates the spatial synchronous detection phase 
measurement in general at a pixel location indicated as 
pixel (x, y,). The video memory 102 is used. The bytes 
20 contained in a window centered at pixel x,y are 

extracted. The window is vertical and may be one pixel 
wide by 25 pixels long. During x-ca 1 i br a t i on , as 
explained in connection with FIG. 5A, the window runs 
horizontally. The bytes in the window are indicated by 
25 block 110. The pixels in the window are convolved with 
sine and cosine attenuated functions which are stored in 
a cosine function memory and a sine function memory 114 
and 116. The functions stored in these memories are 
obtained oy selection of the convolution window width, 
30 the period of sinusoidal activity within the window, and 
a low pass filter function (e.g., Blackman or 
Blackman-llarris) . As indicated above, a 25 pixel width 
is used, however, the width may vary depending upon the 
resolution requirements. The period of the pixels in 
35 the window are then known because of the period of the 
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rjlin3s in the projected grating. For example with the 
grating discussed above, there are generally 2 1/2 to 
3 1/2 periods of the sinewave in a 25 pixel window. 

From this information', sine and cosine 
functions are generated which are pure sine and 
cosines. These functions are multiplied by a digital 
filter function. The Womack article describes the use 
of Hamming filter function. Tt has been found that a 
Blackman filter or Blac kman-Ha r r i s function is 
preferable. The cosine function and low pass filter 
generation processes provide the filtered sine and 
cosine functions ^hich are stored in the memories 114 
and 116 and are convolved with the bytes representing 
the. pixels in the window 110. These are convolutions Ml 
and M2. Ml is the convolution of the pixel information 
shown at 110 which is extracted from the video memory 
102 with the filtered cosine function. The other 
convolution, M2, is with the Blackman f i 1 te r ed. s i ne 
function. The phase is determined by taking the arc 
tangent of tnese convolutions as discussed in the Womack 
article. Ihis phase value is stored and represents the 
phase at. pixel x, y . In this manner, the phase at the 
location of the marker fringe and at the location of the 
window foe pixel 128 is determined. 

Phis absolute phase is converted into 
cumulative phase by measuring the number of 360° 
transitions in phase between the marked fringe location 
and the window location. This is determined by the 
number peaks in the values of the intensity levels along 
the pixel 128 column. Using this cumulative phase 
information and the taole which correlates cummulative 
phase of the marked fringes at each of the 32 planes 
Z Q to 2^ 1# the numbers of the calibration planes 
between which the center of the surface of the object is 
located as well as the distance in the z direction from 
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the origin (2 Q ) is determined, The marked fringe 

table is obtained during calibration by phase 

measurements . upon the calibration target as shown in 

FIG. 15B and as discussed above in connection with 

5 FIGS. 6A and 8. 

As shown in FIG. 15B, the marked fringe is 

located. Then a gray level scaling is applied to the 

pixels in the column along the center line to remove the 

modulation and intensity variation due to the marked 

10 fringe. From the remaining sinusoidal phase 

information, the cumulative phase from the marked fringe 

to the window is measured. This is the phase & 

r mf 

For each of the z planes, the Q c value is stored. 

mt 

These values are used to determine between which 
15 calibration planes (straddling Z planes) the image of 
the surface of the object at the center of the camera 
array 1 ies . 

Next, the phase at all of the pixel positions 
in the window on both sides of the' center pixel column 

20 is calculated using the synchronous detection technique 
discussed above. This provides phase values at each of 
the pixel positions x. This variation in phase for the 
image viewed by camera 3, as shown in FIG. 12A, is shown 
in FIG. 12B. These phase measurements are converted 

25 into absolute distance measurements utilizing the x and 
z calibration tables as explained in connection with 
FIG. 9. 

In the event there are 32 windows of 32 data 
points in each of 32 calibration planes, there will be a 
30 matrix of 32 by 32 by 32 values for each" camera; one 
matrix containing the x distance information and the 
other matrix containing the z phase information. There 
is also a smaller 1 x 32 x 32 matrix containing marked 
fringe information. The x calibration information is 
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obtained using the calibration target as explained in 
general above in connection with FIG. 5. 

FIG. ISA shows how the 32 distance values are 
calculated for each of the windows. First the phase at 
the center pixel is calculated. Then the convolution 
window is moved over in 6 pixel steps and the phase at 

15 additional pixels are calculated. These absolute 
measurements are then converted into unwrapped or 
cumulative phase measurements. The cumulative phase 
measurements are converted to distance since, 1/4 inch 
corresponds to 360° of accummulated phase. The 32 
values ire stored in the memory. The procedure for 
calculating the 15 values start on one side of the 
center pixel (when i=31 starting from 16). Then the 
system is reset back to the center pixel and the 
convolution window is moved in 6 pixel steps until 16 
more x values on the opposite side of the center are 
calculated. There are therefore 32 x distance values, 

16 on one side of the center and 15 on the opposite side 
of the center, plus the center (a total of 32 x 
values). Accordingly, each of the 32 calibration planes 
has 32 oy 32 :< values; thus generating the 32 by 32 by 
32 matricies discussed earlier. 

Similar sets of 32 by 32 phase values are 
scored for each calibration plane in the other aspects. 
There are 32 planes for each camera, each 32 by 32 
matrix may be considered to be a separate look-up table, 
thus for each camera there are 32 x distance look-up 
tables and 32 2 phase look-up tables. 

The 2 phase look-up tables are developed as 
discussed above in connection with FIGS. 7A and 7B. 
FIG. 15C shows how the z phase information at a window 
in a calioration plane is obtained. The procedure is 
modified for the base plane Z=0. At the base plane the 
phase is measured at the center pixel p (pixel 128) 
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and then at each calibration point. It will be noted 
that these calibration points may not necessarily be 
equally spaced at 6 pixels from the center due to the 
bars of the calibration target 1 pattern . If a bar blocks 
5 the calibration point, the phase is measured at a 

unblocked point and interpolated back to the calibration 
point. Then an unwrapped phase difference calculation 
is made to develop the table of 32 values of z phase at 
the window W.. for the base plane. The phase 
10 measurements at each calibration point in the base plane 
Zq is then used for developing the z phase values in 
the next reference plane, % . The temporary memory 
stores the phase values in the preceeding planes so that 
unwrapped phase difference from the previous planes 
15 phase can be determined. The phase values are obtained 
by adding the phase difference in unwrapped phase from 
the previous plane to the phase values in the previous, 
plane for each of the 32 calibration values in the x 
direction. Thus, FIG. 15C shows how the z phase values 
20 are obtained in the Z Q plane as well as the z phase 
values for the same window in each of the other 
calibration planes . 

The total calibration operation is illustrated 
in FIG. 14. For x-ca 1 i br a t ion the calibration target is 
25 illuminated with the constant or uniform illumination 
from the incandescent light sources 46; the projectors 
of the Ronchi pattern being off. Then x-calibration is 
carried out for all windows in the first, Z Q 
calibration plane. Then the projector is strobed and 
30 the marked fringe table is generated for all windows in 
the Z Q plane. During this strobe the video memory has 
sufficient information to obtain the z phase calibration 
for ail windows in the 2 Q plane. Then the stepper 
motor is turned on and the Z 1 plane is calibrated. 
35 Next the stepper motor steps the calibration target to 
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the Z 2 and that plane is calibrated. After all 32 
planes are calibrated the calibration target is rotated 
120° to face the next camera projector assembly and 
reindexea back to the Z Q position. The planes are 
recalibrated for the new camera/projector array. After 
ail 32 calibration planes are calibrated for the second 
camera projector array the calibration target is rotated 
another 120° to face the third camera projector array 
and the 32 planes are again calibrated. All of the 
calioration information is stored in the tables. 

The calibration target can then be removed and 
the system used for object profilometry in three 
dimensions. Once the system is calibrated, it may be 
recalibrated from time to time depending upon changes in 
alignments which may occur during use. For example in 
installations subject to vibration or other changes 
which may affect the positioning of the elements, 
r eca 1 i br a t i on may be needed more frequently than 
otherwise would be the case. 

After the phase is calculated across the window 
from the objects surface these phase measurements are 
converted into x and z distance measurements using the x 
and z phase tables as explained in connection with 
FIG. 9. It will oe noted that the %2 = %? i2 3 wnicn is 
the center pixel in the window. In the event that phase 
measurements are made at pixels other than those spaced 
6 pixels apart, an interpolation program in z and x may 
be used to obtain the calibrated values in the z planes 
for the intermediate pixel positions. The distance 
conversion thus provides a family of x and z values at 
any cross-section of the object under test or at a 
plurality of cross-sections if more than one 
cross-section is to be measured. 

Referring to FIG. 11A, since the object under 
test is imaged onto a number of pixels, there will be a 
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number of x and z values in each aspect. The number of 
X x and 2 1 values equal to the number of pixels in 
camera #1 on which the object is imaged. Similarly 
there will be a number of X 2 and Z 2 values and X 3 
5 and 2 3 values corresponding to the number of pixels in 
the second and third cameras on which the object is 
imaged . 

These values are then compensated for 
temperature variations and also for offsets as shown in 
10 FIG. 11B. The offsets are corrections that are used in 
merging the three cameras* coordinate systems. The 
offsets are measured through the use of a reference 
object, when the measurements are designed to compare 
the measured values with a standard. The standard 
15 object is then measured and the offsets from the known 
standard are stored in an offset table. This standard 
is referred to as a golden master and is measured at a 
reference time and temperature. The offset values are 
indicated as A X Qff and /\ z f f . 
20 T « )e temperature change is also an offset. The ■ 

temperature can be detected by the use of a target which 
remains in the system. This target has a pattern of 
lines similar to that of the calibration target but may 
be smaller, since only a phase measurement at a pixel 
25 position near the top of the CCD sensor array is used. 

This may be pixel 128 in row 12, or closer to the top if 
a smaller convolution window is used. The pixel used 
may be that closest to P 12g in row 12 which is at the 
intersection of bright areas between the projected 
30 grating and the temperature target ruling. Both the x 

and z measurements at this pixel position are used as an 
offset and added to the standard offset in an offset 
accumulator foe each of the three aspects, with the 
offsets added to the values of x and z (or subtracted as 
35 the case may be) to provide new values X^, 2^; 
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X* . , an1 X' , Z' . In order to obtain an 

el — J J 

output in two dimensions, the X' 2 , 2'^ and X'^* 
2' 3 vaij;3s are given +120° and -120° coordinate 
transformations, respectively.* There are therefore a 
succession of x and z values for successive locations 
360° about the cross-section. These x and 2 values are 
stored in memory. 

The aspects overlap each other. In the 
overlapping areas some measurements are inconsistent. 
For example a measurement of greater than 1/8 inch 
corresponds to a phase ambiguity and is considered 
aberrated data. If the data does not change 
moMOt on i c-i i Ly (either increasing or decreasing), the 
probaoi I : !*.y is that the data is serrated. Such data 
points ar.' also removed. It may also be desirable, 
optionally and especially where the centroid of a 
cross-section ls to be computed from the data, to remove 
nil >■ redunoant data point. This will p to vide only data 
in the klms «hicn 00 not overlap. The output data 
rr»pr esent r the c toss-sect i on of the object. It can be 
compared to standard cross-sections, either in area or 
along particular radii tor the purpose of inspection and 
acceptance of parts. The distance measurements may also 
be used to provide machine vision inputs to robot arms 
which mu:;t move precisely with respect to the objects 
under test . 

Siy way of summary the following provides a 
listing of the notation used above and the computations 
carried out by the system in obtaining the absolute 
measurements of dimensional coordinates of the profile 
of an ooject. 

SiZ Calibrated step size. May be any value as 

long as it results in phase changes from Z n 
t d Z n+ j less than 180°. This patent 
application uses £ 2=1/8". 
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Z-CALI B RATI ON 

P Accumulated phase to plane Z n (Just before 

the unknown distance) 

p ' Accumulated phase to plane Z n+1 (Just after 

5 the unknown distance) 

MF-CALI BRATION 

tfttFzn Accumulated phase from Wi to marked fringe 

for plane Z n 

ifMFzn + l Accumulated phase from Wi to marked fringe 
10 for plane 2 n+1 

RUN-TIME MEASUREMENTS 

0123 Phase at center (X=128) for window Wi 

- 0rMF P Accumulated phase from Wi to marked fringe 

&x Phase at pixel, (x # Y-Wi) 

15 U(U X ) Unwrapped phase 
COMPUTATION 

%z Depth between planes Z n and Z n + 1 

%z = u -0 -a ) (i) 

2 Distance 
20 2 = Z n + %2( A z) 

Z Size of the step or incremental distance 

between planes(viz. , aZ = Z-Z'«Z -Z ) 

P 128 Accumulated phase to unknown distance at the 

centre of the video. This is the link between 
25 th e marked fringe calibration and the 

Z-calibration. Since P, P' from the Z 
Calibration are known and %Z was found using 
equation (1), we can find pl28. 

pl28 = P + %z (P'-P) ( 3) 

30 To calculate phase at pixels past the center, 

and then extract the distance, the following 
procedure is performed 
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Calculate Unwrapped phase difference 

U( & Jfr) = bi# x -J? x „ x ) (4) 

Calculate Total phase difference 

Px " Px-1 + U { /i 4)) (5) 

5 Calculate %px 

%Px = (P-Px)/(?-P' ) (6) 

Calculate Distance 

Z x = Z + %p x (Z-Z' ) (7) 

While an exemplary embodiment of the invention 

10 has been disclosed, variations and modification thereof 
within the scope of the invention will undoubtedly 
suggest themselves to those skilled in the art. For 
example, for objects which are not cylindrical the phase 
measurements may be made using temporal synchronous 

15 detection rather than spatial synchronous detection, for 
example translating the projected grating and measuring 
the relative phase shifts at each translation. 
Accordingly trie foregoing description should be taken as 
illustrative and not in a limiting sense. For example, 

20 numerical values may relate to precise values and to 
values substantially the same as the precise values* 
The present application includesthe disclosures of 
the accompanying claims, drawings, and abstract of 
the present invention, and those disclosures are 

25 hereby incorporated in to the present description. 
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CLAIMS 

1, A system or apparatus for measuring the profile 

of an object by digital moire interf erometry which 
comprises means for calibrating' a volume in space in 
5 which volume the object is disposed during measurement 
to provide reference phase information identifying the 
locations of a plurality of sets of positions in said 
volume, each of said sets being in a different one of 
a plurality of parallel calibration planes spaced from 

10 each other in said volume; means for obtaining phase 
information from the profile of said object in said 
volume at said plurality of positions; means for 
identifying from the phase information obtained from 
said profile the ones of said calibration planes on 

15 opposite sides of said profile; and means for 

determining from the differences between said phase 
information from said profile and said sets of phase 
information for said calibration planes on said 
opposite sides of said profile absolute distance 

20 measurements of the dimensions of said profile, 
at least through a cross-section thereof. 
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2. The system or apparatus according to claim 1 
wherein said object is three dimensional and said 
calibrating means includes means for calibrating said 
volume in a plurality of angularly offset aspects of 

5 said volume, a different plurality of said parallel 
planes in each of said aspects being perpendicular to 
radial lines through an axis of rotation, the angular 
offset being the angle between said radial lines. 

3. The system or apparatus according to claim 1 
10 wherein said calibrating means comprises a target having 

a pattern of lines spaced from each other and 
paralleling a first of three orthogonal axes; means 
for illuminating said target, means for deriving a 
plurality of first tables each containing data 

15 representing a matrix of distance values derived from 

the phase variations in the intensity of light reflected 
from said target at said set of positions from each of 
said plurality of said planes, said set of positions 
being along a second of said orthogonal axes which is 

20 also perpendicular to said lines of said target, said 
planes being spaced along a third of said orthogonal 
axes, each of said tables being derived with said 
target pattern disposed in a different one of said 
plurality of parallel planes; means for projecting a 

25 ruling having a plurality of spaced parallel lines 
upon said target pattern so that the lines of said 
ruling are perpendicular to the lines of said target; 
means for optically marking at least a predetermined one 
of said ruling lines; means for deriving a table 

30 containing data representing the difference in the 
phase in the intensity of the light reflected from 
said marked line and at a predetermined position in a 
window spaced at a predetermined position perpendicualr 
to the pro jected lines, said last named table 

35 identifying the positions of each of said planes along 
said third of said orthogonal axes; 
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means for deriving a second plurality of tables each 
containing data representing a matrix of phase 
information from light reflected from said target 
pattern with said ruling projected thereon at said set 
5 of positions in each of said plurality of planes; 

and said means for determing said dimensions comprising 
means for deriving said dimensions of said profile along 
said second axis from those of said first plurality of 
tables for planes straddling said profile and from the 
10 phase information from said object, and means for 

deriving said dimensions along said third of said axes 
from those of said second plurality of tables for the 
planes straddling said profile and from the phase 
information from said object. 
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4- . The system or apparatus according to 
claim 3 wherein said calibrating means and said means 
for obtaining phase information from said profile of 
said object comprises a camera having a sensor defining 

5 a matrix of pixel positions along said first and second 
axes, means for deriving digital data signals from each 
of said positions, means for convolving said digital 
data from those of said positions in successive areas, 
each constituting a convolution window, with digital 
10 signals representing sine and cosine functions to 

provide digital outputs, and means for deriving said 
phase information in accordance with said digital 
outputs. 

5- The system or apparatus according to 
15 claim 4 wherein said optical marking means includes 

means for projecting a marked area at said predetermined 
one line. 

6. The system or apparatus according to 

claim 5 wherein said marking means includes means for 

20 projecting a plurality of said marked areas, one at said 
predetermined line and another at a line spaced a 
predetermined number of lines from said predetermined 
one line, each of said areas having different intensity 
of i 1 lumi nation. 

25 7» The system or apparatus according to 

claim 5 wherein said calibrating means and said means 
for obtaining said phase information the profile of said 
object comprises a camera definining a matrix of pixel 
positions along said first and second axis, means for 

30 deriving digital data signals from said positions, means 
for processing said digital data signals to remove the 
variation in value thereof at positions corresponding to 
said marked area, and means for obtaining said phase 
information from said digital data signals from which 
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said variation in value due to said marked area has been 
removed . 

8. The system or apparatus according to 

claim 3 wherein said calibrating means includes means 
5 for deriving said phase information in terms of 

unwrapped phase having means for accumulating phase at a 
preceding position, 360° of phase times the number of 
lines between said preceding and a succeeding position, 
and. the phase at said succeeding position. 

10 9* Tne system or apparatus according to 

claim 3 further comprising means for deriving offset 
distance measurements in said volume, and means for 
correcting the distances derived from said profile of 
said object in accordance with said offset distance 

15 measurements. 

10. The system according to claim 9 wherein 
said offset distance measurements deriving means 
includes another target in said volume adjacent the 
object under test, said other target having a pattern of 

20 lines thereon, means for illuminating said other target 
in said volume to obtain offsets corresponding to 
temperature variation, and whereby said offset distance 
deriving means includes an object of reference profile, 
and the means set forth in claim 3 for deriving said 

25 distances. 

11. A system or apparatus for measuring the 
profile of an object by digital moire inter ferometry / 
substantially as hereinbefore described with reference 
to the accompanying drawings. 

30 12. A method of digital moire prof i lomet ry , 

comprising utilizing at least one system or apparatus 
according to any one of the preceding claims. 
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Amendments to the claims 
have been filed as follows 



13. A method for measuring the 

profile of an object by digital moire i nter f e r omet r y 
which comprises tiie following steps: calibrating a 
volume in space- in which volume the object is disposed 
5 during measurement to provide reference phase 

information identifying the locations of a plurality of 
sets of positions in said volume, each of said sets 
being in a different one of a plurality of parallel 
calibration planes spaced from each other in said 

10 volume; obtaining phase information from the profile of 
sal: juivci. in saii volume at said plurality of 
po.; : l i »i:i?5 ; i-s _*.!:_ 1 1 / 1 ng from the phase information 
obtained trom'sdi j profile the ones of said calibration 
planes on opposite sides of said profile; and 

15 determining from the differences between said phase 
information from said profile and said sets of phase 
information for said calibration planes on said opposite 
sides of said profile absolute distance measurements of 
the dimensions of said profile, at least through a 

20 cross-section thereof. 

14. The method according to 

claim 13 wherein said object is three dimensional and 
said c j : i -;. »; -J: i ng step is by repeating said calibrating 
step for a plurality of angularly offset aspects of said 
25 volume, a different plurality of said parallel planes in 
each of said aspects being perpendicular to radial lines 
through an axis of rotation, the angular offset being 
the angle between said radial lines. 

15. The method according to 

30 claim. 13 wherein said calibrating step is carried out by 
illuminating a target having a pattern of lines spaced 
from each other and paralleling a first of three 
orthogonal axes; deriving a plurality of first tables 
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each containing data representing a matrix of distance 
values derived from phase variations in the intensity of 
light reflected from said target at said set of 
positions from each of said plurality of said planes, 
said set of positions being along a second of said 
orthogonal axes which is also perpendicular to said 
lines of said target, said planes being spaced along a 
third of said orthogonal axes, each of said tables being 
derived with said target pattern disposed in a different 
one of said plurality of parallel planes; projecting a 
ruling having a plurality of spaced parallel lines upon 
said target pattern so that the lines of said ruling are 
perpendicular to the lines of said target; optically 
marking at least a predetermined one of said ruling 
lines; means for deriving a table containing data 
representing the difference in the phase in the 
intensity of the light reflected from said marked line 
and at a predetermined position in a window spaced at a 
predetermined position perpendicular to . the projected 
lines, said last named table identifying the positions 
of each of said planes along said third of said 
orthogonal axes; deriving a second plurality of tables 
each containing data representing a matrix of phase 
information from light reflected from said target 
pattern with said ruling projected thereon at said set 
of positions in each of said plurality of planes; and 
said step of determining said dimensions comprising the 
steps of deriving said dimensions of said profile along 

said second axis from those of said first plurality of 

> 

tables for planes straddling said profile and from the 
phase information from said object, and deriving said 
dimensions along said third of said axes from those of 
said second plurality of tables for the planes 
straddling said profile and from the phase information ' 
from said object. 
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